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Abstract 

We present an updated theoretical prediction of e'/e, using the next-to-leading AS — 1 effective 
hamiltonian and lattice QCD matrix elements. The CP violating phase is costrained by using both 
the experimental values of e and Xd, assuming the theoretical determination of fg. Predictions of 
cos 5 and sin 2/3 are also obtained in this way. For e'/e, our estimate is e'/e = (2.8 ± 2.4) x 10~ 4 . 



We repeat the combined analysis of the CP violation 
parameter e and the B-mixing parameter Xd in order to 
estimate e'/e, along the lines followed in refs. JjL |^| . The 
main steps of this analysis are the following: 

1 The CP violating phase 5 of the CKM matrix is 
constrained by comparing the theoretical prediction 
for e with its experimental value. To this purpose, 
the relevant formula is 

|e|c=o = C e B K A 2 \ 6 asmS{F(x c ,x t ) + (1) 
F(x t )[A 2 X A (l-acos8)]-F(x c )}, 

where x q — m 2 /M^r and the functions F(xi) and 
F(xi,Xj) are the so-called Inami-Lim functions ||, 
obtained from the calculation of the basic box- 
diagram and including QCD corrections. F(xt) is 
known at the next-to-leading order, which has been 
included in our calculation W&. In eq. (hi), 
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where AMr- is the mass difference between the 
two neutral kaon mass eigenstates. Moreover, p = 
a cos 8 and rj = a sin 6, where A, A, p and rj 
are the parameters of the CKM matrix in the 
Wolfenstein parametrization H. Finally, Bk is the 
renormalization group invariant B- factor |J . 
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The theoretical estimate of the B-meson coupling 
constant fs is used to further constrain S, by 
comparing the theoretical prediction of the B mixing 
parameter Xd with its experimental value. From the 
AB — 2 effective hamiltonian, one can derive 
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where Bb is the B-parameter relevant for B — B 
matrix element 



(B d \(dlib) 2 \B d ) = 8/3f 2 M 2 B B . 
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Notice that /bB b must be known, for the 
experimental value of Xd to give a constraint on S. 

3 From the AS = 1 effective hamiltonian, one 
can calculate the expression of e' in terms of 
CKM matrix elements, Wilson coefficients and local 
operator matrix elements. One has 
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where (Im^)' and ItjuAq are given by 
W = -G F Im(v*V td ){-(C 6 B 6 + 
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and 



(W 2 )' = -G F Im(y t * s V td ) \c 7 B* /2 (I- 



6 



(7) 



C 9 B 3 9 /2 =^ 



,2X 
"3 



Here a; = Rey^/Re/lo and we have introduced 
(Im^)' defined as 



lmA 2 = (ImA 2 y + Oj B (wImAo). 



(8) 



fi/s accounts for the isospin breaking contribution, 
see for example ref. [Q. The Wilson coefficients Ci 
have been evaluated at the next-to-leading order for 
/i = 2 GeV, using the anomalous dimension matrices 
£] and the initial conditions 



given in refs. 
computed in refs. [[L0, ll[ (given for HV in ref. ]12fl). 
Concerning the local operator matrix elements, their 
values are given by a set {Bi} of B-parameters 
multiplied by the vacuum insertion approximation 
results. In turn, these can be written in terms of the 
three quantities (see eq. (||) and eq. (^)) 
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The numerical values of the £?-parameters have been 
taken from lattice calculations For those B- 
factors which have not yet been computed on the 
lattice, we have used educated guesses, see ref. Jl3[ . 

More details on this analysis can be found in ref. 
H 0- Compared to our previous work 0, the main 
improvements are the following: 



Parameters 


(fBB' L J 2 ) th = (200 ±40) MeV 
m s (2GeV) = (128 ± 18) MeV 

A QCD = ( 330 ± 10 °) MeV 
t b = (1.49 ± 0.12) x 10~ 12 s 
\V ub /V cb \ = Act = 0.080 ± 0.015 


V cb = A\ 2 = 0.040 ± 0.006 
x d = 0.685 ± 0.076 
n IB = 0.25 ± 0.10 
m t = (174 ± 17) GeV 



Table 1. Values of the fluctuating parameters used in the 
numerical analysis. 



The constraint on S coming from x d is used in the 
analysis, taking Jb from the theory. Since there 
is increasing theoretical evidence that the value of 
/b is large (~ 200 MeV) and that the relevant B- 
parameter Bb is close to one, this constraint is quite 
effective. 

Updated values of the experimental parameters 
entering in the phenomenological analysis, such 



as the B meson lifetime tb, the B^-B^ mixing 



parameter x d , the CKM matrix elements, (|Vd,|, 
|Kifc|/|Kb|), etc., have been used. 

3 The value of the strange quark mass m s has been 
taken from lattice calculations jn|, thus making 
a more consistent use of lattice results for the B- 
parameters of the relevant penguin operators. 

4 All the results are presented with an estimate of the 
corresponding errors. These errors come from the 
limited precision of measured quantities, e.g. tb, 
and from theoretical uncertainties, e.g. the values of 
hadronic matrix elements. 

The results of our analysis have been obtained by 
varying the experimental quantities, e.g. the value of the 
top mass m t , tb, etc. and the theoretical parameters, 
e.g. the B-parameters, the strange quark mass m s (fi), 
etc., according to their errors. Values and errors of the 
input quantities used in the following are reported in 
tables |-|. We assume a gaussian distribution for the 
experimental quantities and a flat distribution (with 
a width of 2a) for the theoretical ones. The only 
exception is m s (/i), taken from quenched lattice QCD 
calculations, for which we have assumed a gaussian 
distribution, according to the results of ref. fl6|| . 

The theoretical predictions (cos S, sin 2/3, e'/e, etc.) 
depend on several fluctuating parameters. We have 
obtained numerically their distributions, from which 
we have calculated the central values and the errors 
reported below. 

Using the values given in the tables and the formulae 
given previously, we have obtained the following results: 

1 The distribution for cos 5, obtained by comparing 
the experimental value of e to its theoretical 
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Figure 1. Distributions of values for cos<5, sin 2/3 and e'/e, using the values of the parameters given in the tables. 
The solid istograms are obtained without using the constraint. The dashed ones use this constraint, assuming 
1/2 

that fsBg = 200 ± 40 MeV. The contour-plot of the event distribution in the p — rj plane is also shown. 



Constants 


G F = 1.16634 x 10~ b GeV-^ 


f„ = 132 MeV 


m c = 1.5 GeV 


} K = 160 MeV 


mt = 4.5 GeV 


A = sin6» c = 0.221 


M w = 80.6 GeV 


e exp = 2.268 X 10" 3 


Mir = 140 MeV 


ReA = 2.7 x 10" 7 GeV 


M K = 490 MeV 


o = 0.045 


M B = 5.278 GeV 


p = 2 GeV 


AM K = 3.521 x 10~ 12 MeV 





Table 2. Values of the constants used in the numerical analysis. 

prediction, is given in figure [l]. As already noticed 
in refs. [jl], |) and |lj], [l8|], large values of fg and m t 
favour cos 8 > 0, given the current measurement of 
x d . When the condition 160 MeV < f B B^ 2 < 240 
MeV is imposed (/^-cut), most of the negative 
solutions disappear, giving the dashed istogram of 
figure [I], from which we estimate 

cos 8 = 0.47 ± 0.32 . (10) 

2 The value of sin 2/3 depends on cos<5. The 
distribution of sin 2/3 is shown in figure [l], without 
(solid) and with (dashed) the /s-cut. When the Jb~ 
cut is imposed, one gets larger values of sin 2/3 [Q. 



From the dashed distribution, we obtain 

sin2/3 = 0.65 ± 0.12 . (11) 

Figure |l| also contains the contour-plot of the event 
distribution in the p — rj plane, showing the effect 
of the e and Xd constraints, when the /s-cut is 
imposed. 

3 In figure ||, several informations on e'/ e are provided. 
Contour-plots of the distribution of the generated 
events in the e'/e-cos 8 plane are shown, without 
and with the /s-cut. One notices a very mild 
dependence of e'/e on cos 8. As a consequence one 
obtains approximatively the same prediction in the 
two cases (see also figure [j]) 

e'/e = (2.3 ± 2.1) x 10~ 4 no - cut, (12) 

and 

e'/e = (2.8 ± 2.4) x 10~ 4 f B - cut. (13) 

In figure^, we also give e'/e as a function of mt- The 
band corresponds to the 2a prediction. 

In spite of several differences, the bulk of our results 
overlap with those of ref. It is reassuring 

that theoretical predictions, obtained by using different 
approaches to evaluate the operator matrix elements, 
are in good agreement. 
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Figure 2. Above, contour-plots of the event distributions in the plane e'/e-cosS without and with the /s-cut. 
Below, e'/e as a function of mt- 



B-parameters 


B K = 0.75 ±0.15 
B{_. 2 = 0-0.15<» 
#5,6 = Bf_ll ] =1-0 ±0.2 


Bf /2) = 0.62 ±0.10 
B- iA = 1 - 6 W 

-°7-8-9 1 



Table 3. Values of the B-parameters, for operators 
renormalized at the scale fi = 2 GeV. The only exception is Bjf, 

3/2 

which is the RG invariant B-parameter. B g has been taken 
equal to Bk, at any scale. The value reported in the table is 
Bg /2 ([i = 2GeV). Entries with a W are educated guesses, the 
others are taken from lattice QCD calculations. 

On the basis of the latest analyses, it seems very 
difficult for e'/e to be larger than 10 x 10~ 4 . This may 
happen by taking the matrix elements of the dominant 
operators, Qe and Qg, much different than usually 
assumed. One possibility, discussed in ref. E3), is to 
take B 6 ~ 2 and B s ~ 1, instead of the usual values 

~ Bg ~ 1. To our knowledge, no coherent theoretical 
approach can accomodate so large value of Bq. 
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